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EQUILIBRIUM  OF  A  TOROIDAL  PINCH  IN  A  MAGNETIC  FIELD 


V.  D.  Shafranov 


The  conditions  of  equilibrium  of  an  axisymmetric 
toroidal  pinch  are  obtained  by  analysis  with  respect 
to  a/R  (a  and  R  —  the  small  and  large  radius  of  torus) 
without  concrete  representation  of  the  type  of  current 
distribution  and  of  longitudinal  magnetic  field  with 
respect  to  the  pinch  profile.  It  is  assumed  only, 
that  in  a  zero  approximation  (R  =  «)  the  pinch  has  a 
cylindrical  symmetry. 

Refined  formulas  are  given  for  pinch  displacement 
In  a  conducting  housing  that  are  suitable  in  the  second 
approximation,  i.e.,  taking  into  account  terms  of  the 
order  b2/R2 ,  where  b  i3  the  housing  radius. 


Introduction 

The  equilibrium  conditions  of  a  toroidal  axisymmetric  pinch  in  a 
magnetic  -field  were  derived  in  earlier  articles  [1-6]  under  certain 

i 

hypothesis  of  the  current  distribution  with  respect  to  a  pinch  cross 
section.  In  previous  articles  [1-5]  idealized  configurations  either 
with  a  surface  current  or  with  a  uniform  current  were  examined.  One 

article  [6]  cited  the  solution  for  the  case  where  the  surface  current 

i  ? 

distribution  with  respect  to  the  cross  section  was  described  as  a 

:  ,  ' 

zero  Bessel  function  (from  an  actual  or  imaginary  orgu-vent) .  .  Below, 

•  •  f.  > ; C. ,■;<>£ 
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the  equilibrium  conditions  are  obtained  without  concrete  representa¬ 
tion  of  the  type  of  current  density  distribution  or  longitudinal  mag¬ 
netic  field  with  respect  to  the  section.  Presented  only  is  that 
a  «  R  and  that  in  the  zero  approximation  (R  =  «)  the  distribution  la 
cylindrlcally  symmetrical. 

To  find  the  magnetic  field  confining  the  pinch  in  equilibrium, 
we  must  know  the  field  distribution  on  the  pinch  surface.  Thus,  for 
example,  in  an  idealized  case  of  a  surface  current  and  with  no  lon¬ 
gitudinal  magnetic  field,  the  pressure  of  the  plasma  in  the  pinch  is 
constant  p  =  const,  and  consequently  the  meridional  magnetic  field 
at  the  pinch  surface  is  constant.  This  condition  permits  us  to  find 
the  general  magnetic  field  outside  of  the  plasma  [2],  In  the  case  of 
space  currents,  the  meridional  magnetic  field  at  the  pinch  surface  is 
already  variable,  but  when  «  1  it  can,  apparently,  be  expressed  in 
terms  of  such  characteristics  as  internal  self-induction  of  a  unit  of 
pinch  length  and  average  pressure  p  of  the  plasma  with  respect  to 
the  section.  For  the  surface  of  the  pinch  with  distributed  current 
we  can  assume  an  arbitrary  magnetic  surface  on  which  the  current  den¬ 
sity  is  neglibible  as  compared  with  the  average  current  density  in 
the  pinch.  In  the  zero  approximation  {R  =  «)  magnetic  surfaces  are 
a  system  of  concentric  cylinders  (p  =  const) .  In  the  first  approxima¬ 
tion  with  respect  to  a/R  the  magnetic  surfaces  are  a  system  of  enclosed 
toruses  whose  cross  sections  remain  as  circles  but  are  no  longer  con¬ 
centric.  • 

In  this  article  the  magnetic  field  distribution  over  magnetic 
surface  is  derived  in  part  2  by  examining  the  balance  of  forces  acting 

on  a  volume  element  of  wedge-shaped  form,  formed  by  the  rotation  of  ! 

1  ' 

a  ..ector  of  the  section  of  magnetic  surface  around  the  axis  of  1  ' 

-  i  .  o  ..  j 
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symmetry  of  the  torus.  The  boundary  condition  (21)  derived  In  that 
section  Is  used  In  part  5  for  calculation  of  the  external  magnetic 
field  needed  to  confine  an  annular  pinch  in  equilibrium.  In  part  4, 
using  the  same  condition  (21) ,  we  obtain  the  formulas  describing  the 
equilibrium  of  a  toroidal  pinch  inside  a  conducting  housing.  The 
formulas  for  magnetic  field  distribution  outside  the  plasma  and  for 
the  displacement  value  of  the  center  of  the  plasma  section  relative 
to  a  section  of  Die  housing  are  presented.  We  will  examine  the  dlstor 
tion  of  the  shape  of  its  pinch  section  associated  with  the  toroidal  ' 
effect  and  the  influence  of  various  factors  (extraneous  factors, 
”scatteredw  fields,  separations  in  the  housing,  the  finite  electrical 
conductivity  of  the  latter)  on  the  equilibrium  of  the  pinch. 

2.  Magnetic  Field  Distribution  on  a  Magnetic  Surface 


The  equation  of  equilibrium  of  plasma  in  a  magnetic  field 


(i) 


can  be  written,  a3  i3  known,  by  means  of  the  tensor  of  maxwellian 
stresses  in  an  integral  form  1 

=  0, 

where 


r«fi=(p+Wj* 


4* 


(3) 


I 

s  •  >  - 

Here  integration  in  (2)  13  performed  with  respect  to  the  surface 

} 

;•  -.‘bitrarily  bounding  an  isolated  volume.  To  determine  the  magnetic 

I 

field  on  a  magnetic  surface,  whose  section  radius  equals  a,  we  apply 

•  I  ' 

relation  (2)  to  the  volume  element  which  is  bounded  by  coordinate 

’(  r  .*  \  *  #*  }*  '.rf: 
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surfaces  (Fig.  1) 


9  -t-  dip  =a  const;  0  <  Q  <  a. 


(*> 


The  coordinates  p,  9,  <u  are  connected  with  the  cylindrical  coordinate# 
r,  9,  z  by  the  relation 


r  —  R  +  Q  cos  (1);  z  =  Qsin<i>. 


(5) 


where  R  is  the  distance  of  the  center  of  the  section  of  the  examined 
magnetic  surface  from  the  axis  of  symmetry.  The  corresponding  unit 
vectors  are  determined  by  the  relations 


eQ  —  cr  cos  (o  t-  e,  sin  <*>, 
c<#  =  e,  cos  a)  —  er  sin  co. 


f**  =  _e 

du  e° 


— 


(6) 


We  will  carry  out  the  calculation  in  the  first  approximation  of 
expansion  with  respect  to  a/R.  In  the  zero  approximation  (R  =  »)  we 
will  deal  with  a  cylindrical  pinch,  for  which  the  relation  [71 


-  ,  bI  nl 


(7) 


is  valid. 

Here  B.^  and  Bg  are  the  internal  and  external  longitudinal  fields; 


B.  = 


2J 


13  the  current  field:  pa  is  the  pressure  where  p  =  a.  The  line 

cl  Ca  a  r 

o 

indicates  averaging  over  the  pinch  section  0  <  p  <  a.  In  the  first 
approximation  terms  of  the  type  b°b(1)  appear  in  tensor  T^,  wherein 
B^1)  ~  ®.  Since  on  averaging  over  a  section  of  the  toroidal 

pinch  such. terms  disappear,  then  (7)  holds  true  al30  in  the  first 
approximation  if  B  is  taken  a3  the  value  of  the  longitudinal  magnetic 
field  on  a  magnetic  surface  with  radius  a  for  oj  =  ~.  As  is  known,  in 

j  “  «  • 

an  axlsynudotrical  pinch  the  longitudinal  magnetic  field  changes  along. 

f  r  .  'h  .  .  .  S  ft  r ’  .  1  - 
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the  magnetic  surface  i>  p,  o  =  const  with  respect  to  the  law  [1] 

„  2/  (\|>)  2/  (>w  { *  q  ,  \ 

B*  -  -  ~tjr  C 1  ~  T  cos  «  +  •  •  0  - 

•  =  #.((?)(  1  — -jpcosw-f  . 


(8) 


The  pressure  of  plasma  on  the  magnetic  surface  maintains  a  constant 
value  p  =  p  M. 

We  will  designate  by  dS  ,  dS  ,  dS  the  surface  elements  the  nor-. 

r  * 

mals  to  which  are  unit  vectors  e  ,  e  .  The  resultant  of  forces 

P '  9  CD 

acting  on  the  surfaces  dS  (9)  and  dS  (9  +  d9)  =  dS  (9)  =  p  dp  do,  is  ! 

'  '  »  j 


•  de 

F**  =  -  d«pd(D^  T^a-j^QdQ, 


(9) 


where  summation  with  respect  to  a  =  p,  9*  o  is  Inferred.  This  is  a 

force  of  the  first  order  of  smallness,  therefore  tensor  T  •  should  be 

90 

calculated. only  in  the  zero  approximation: 


!•  ..r  —  FUJ  =*  —  d< p  du>  ^  (e,.  sin  »  — 

-  oQ  cos  to)  —  T * uc9  sin  o>]  q  dp. 


(10) 


The  resultant  of  forces  acting  on  surfaces  dS^fo)  and  dS^fco .  +  dco)  is 


analogously  determined.  Here  we  need  only  calculate  that  on  changing 

'  '  ' .  i  i 

from  o  to  o  +  do  both  the  unit  vectors  and  T^  and  al3o  dS^  =  (R  + 

+  p  co3  o)  dp  d9  change.  We  will  present  the  result: 


a 

F..=  F*.  +  (  /?e#  ^  7^  dp)  + 

0 

+  rfo)  |  c0  ^  COS  (I)  ^  dQ  + 

-s- « S  (rsi- 

-r  ^  sin  (o  [TU,QdQ-n  J  (nj>+ 

+  +•*  [sin<0  \  TttfdQ- 

-«$#*>]}• 


J  .. 


uc?  i-i*e 


'  (11)  i 
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Finally,  taking  into  account  that  B  (a)  =  0,  and  consequently  T  (a) 

p  pm 

=  Tp^(a)  =  0,  we  will  find 


F<o  —  —  eQT qq(R  -:,  a  cos  w)  a  d<f  c/cj i 
==««[  —  4-  4- 

+  r>*cosu)]3FV+Fii'. 


(12) 


Setting  the  sura  of  these  three  forces  (10)-(12)  equal  to  zero,  we  will 
obtain  four  conditions  of  equilibrium  of  an  isolated  volume  element: 
equilibrium  with  respect  to  p  in  the  zero  approximation  and  equilibrium 
with  respect  to  9,  cd  and  p  in  the  first  approximation: 


^  Tlu dQ—Tlj (a) a;. 

(13) 

~~[T‘"dQ=2-^yn<QdQ-, 

(11) 

T  -jj*-)  dQ  -  \  (C„  -  TU)  Q  dQ-, 

(15) 

COS  (0  \ 

\  (TU  +  'flu)  QdQ-R  $  -  7™  J  X 

(16) 

X  cIq  —  Tl<p}  cos  <0  —  7 '$Ra  **  0. 

The  equilibrium  condition  in  the  zero  approximation  (lj5)  can  be  obtained 

after  having  averaged  (l)  over  the  plasma  section.  Conditions  (14), 

/  X  d2rp(J  ' 

(15)  can  be  written  differently  if  we  consider  that  — g^z- 


aaT(l>  =  _T(x)# 


since  ~  cos  aa  *  Thu3>  ^or  example,  (15)  can  be  written  as 

5  (^-  -  r”0 dQ  $ (TL  -  QdQ- 

Eliminating  the  second  integral  in  (16)  with  the  help  of  (17)  we  will 

'f w  (a)  =■  -*f  fos «  [  —  (<*)  + 

« 

+  ~ r  |  dQ  J 


(18) 


obtain 


In  the  examined  approximation 


7## 


rw 


wi  _i£ 
&T“  8^ 


(19) 


(the  term  B*/8ir  is  dropped  3ince  it  is  of  the  second  order  of  small¬ 
ness) .  Using  (8) ,  we  will  find 


rtf,’ (*) 


DaBX _ a? 

4n  8ft 


2 


-y*  COS  10. 


(20) 


From  (l8)-(20)  and  (?)  we  will  obtain  the  unknown  meridional  field  on 
the  magnetic  surface  with  radius  a:  B^a)  =  +  B^*^ ,  where 


B 


-®-Bb  coso> 


r  it  .  to  (p—pj 

L  2  +  Bl 


Here  the  internal  self-induction  of  a  unit  plasma  length  is 


(21) 


1^2)  ; 

I 

(for  surface  current,  for  example,  it  is  l ^  =  0;  for  uniform  current 
it  is  =  1/2).  | 

We  note  that  for  the  derivation  of  (18)  a  concrete  representation 
of  tensor  To  was  not  used;  only  the  condition  T  (a)  =  T  (a)  =  0, 

up  r  * 

was  used,  therefore.  Expression  (18)  is  applicable  not  only  to 

1  i 

hydromagnetic  equilibrium  configurations.  If,  for  example,  we  take 
Ta0  =  pSa£  +  as  the  tensor  then  Condition  (18)  permits  us 

to  find  a  hydrodynamic  flow  originating  about  an  annular  vortex.  ' 

Relation  (21)  is  correct  for  any  magnetic  surface  In  an  equilib¬ 
rium  configuration,  whose  section  radius  is  a  «  R.  Subsequently  ! 
this  condition  will  be  applied  to  anpineh  surface, n  i.e.,  applied  to 
that  magnetic  surface  outside  of  which  the  magnetic  current  density  • 

’ :  1  >  J 

can  be  taken  as  equal  to  zero.  t  •  ' 
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5.  Equilibrium  Without  External  Conductors 


A  meridional  magnetic  field  outside  of  the  pinch  can  be  written 


as 


B— 


2.V 


(23) 


wherein  the  lines  of  magnetic  force  are  arranged  on  the  surface 
^  =  const.  When  a  «  R  [4] 


where 


lj/  ■  [2(eosh  <►—  cos<l))]~1,*  X 
X[/’oW+2f,  (0)  cos®]. 


Fo  (0) = aoS»  (0)  +  V#  (♦); 


(24) 


Here  a0,  b0,  aXj  bi  are  the  constants  subject  to  determination;  g0ji 
gi,  fo>  fi  are  the  Fok  toroidal  functions  [8];  <o  are  the  toroidal 
coordinates  which  when  p  «  R  are  connected  with  p  and  co  by  the 
relations 


4’  e~*  =  1r +-&('-&')  cosu>- 


(25) 


When  p  «  R 


\J  -  b#  (d -f-  2  Id  2 — 2)  + 

+  {[«.+  4i.^  +  21n2-2)]  «"*  + 

+  tje*— cos®. 


,  re 


The  confining  field  ^  is  obtained  when  bo  =  bi  »  0: 


(26) 


4'”“^+  (fl0  — ai)c”*  C0SG)s 


(27) 


;rr.'  '<r- 
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The  corresponding  magnetic  field  la  obtained  by  the  formulas  i 

i/*  n.  k r*  1  *N'r>  wr  (28] 

=0,  Bt  -  dr  inli*  *  '<0'  , 

The  constants  a0,  bo,  aj,  bi  determined  from  the  conditions  of  con-  i 
stancy  and  continuity  on  the  pinch  surface,  are  (assuming  f  m  ^ 
=  0)  equal  to 

•  *  i 

);  (29J. 

.  3«  «* 

®‘  =  *i6^a«-  i 

Thus,  In  agreement  with  Expressions  (28)  and  (29),  the  external  mag-  ' 

netlc  field  needed  to  confine  the  plasma  In  equilibrium  Is  1 


(50): 

T 

In  particular,  from  this  formula  the  values  of  B^e  found  In  earlier 

z 

articles  [1-6]  are  obtained.  Assuming  B^e  =  0,  from  (30)  and  (7)  we 

z 

will  find  the  conditions  of  equilibrium  of  a  magnetic  field  in  a 
toroidal  space  in  the  case  where  the  magnetic  field  Immersed  in  plasma 
is  analogous  to  the  magnetic  field  in  sun  spots  • 

,  1 

8ji/>,  ~  +  /n  ( In  H  _  i.  V  1 

(31)  | 

'.*»-*(  +  '! 

i 

In  a  limiting  case  of  surface  current  (l^  =  0,  p  =  const)  these  formulas 

agree  with  the  results  obtained  in  a  previous  article  [1).  | 

•  ,  !  .  5  '  . ; 

4.  Equilibrium  In3lde  a  Conducting  Housing  ( 

'•  ■  ,  l 

Important  In  practice  Is  the  situation  of  equilibrium  in  a  con-  ] 

l  *  .  2  ■  -  . 

ducting  housing  which  13  used,  for  example,  in  an  experimental "Tokamak" 

'  ••  •  .  i  '  •  ' 

I  >  .  ■ . I  . - . .  ..  , . 1  ..  .  i 
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system  Intended  to  create  and  study  high-temperature  plasma.  The  forces 
resulting  from  the  removal  of  the  line  current  from  the  housing  center 
wa3  calculated  by  M.  A.  Lentovich  [9,  10].  In  one  of  the  articles  [9J 
cases  are  examined  where  the  thickness  of  skin  layer  is  significantly 
less  and  significantly  greater  than  the  housing  wall  thickness.  Another 
article  [10]  took  into  account  influence  of  the  housing  Joints  needed 
to  lead  the  electrical  field  into  the  discharge  chamber.  These  articles 
[9,  10]  solve  the  principal  problem  of  the  possibility  of  confining 
a  pinch  inside  a  toroidal  conducting  housing.  By  using  the  expression 
for  forces  obtained  in  these  articles  [9*  10],  we  can  estimate  the  dis¬ 
placement  value  A  of  the  pinch  inside  the  housing  [11].  The  experi- 

1 

ments  on  "Tokamak"  showed  that  certain  substantial  peculiarities 
of  discharge  are  apparently  associated  with  a  shift  of  the  pinch  during 
discharge  which  leads  to  contact  of  the  plasma  with  the  diaphragm 
confining  the  discharge  [12].  Thus  the  problem  arises  of  finding  a 
more  precise  value  of  the  pinch  displacement  in  the  housing.  Expres¬ 
sion  (21),  which  was  found  above  for  the  meridional  field  on  the  pinch 
surface,  yields  the  solution  to  thi3  problem.  Below,  expressions  are 

i 

obtained  for  pinch  displacement  in  a  conducting  housing  in  the  second j 
approximation  of  expansion  with  respect  to  b/R  and  also  determined 
is  the  distortion  of  the  shape  of  the  pinch  section  caused  by  addi¬ 
tional  so-called  scattered  fields.* 

In  the  second  approximation  of  expansion  with  respect  to  1/R,  the 

i 

pinch  section  will  no  longer  be  round  but  will  have  an  elliptical  form 

,  •  1 

I 

Q  —  a  +  6  cos  2&). 

3  ! 

-  I 

1 

*  L.  A.  Artsimovich  acquainted  the  author  with  the  last  problem* 

1  1 '  . :  1 . .  •  .  0 .  .  .  ; 
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We  will  assume  that  the  housing  section  in  the  general  ca3e  also  has 
an  elliptical  form.  #  In  a  polar  system  of  coordinates  with  an  origin 
in  the  center  of  the  housing  section,  let  the  equation  of  the  surface 
of  the  latter  be 


q'  a  b  -f  x  cos  2«'. 


(33) 


Then  in  the  system  of  coordinates  p,  to  related  to  the  pinch,  whose 
cross  section  center  is  shifted  from  the  center  of  the  housing  section 
(outward)  to  a  distance  A,  the  equation  of  the  housing  surface  has  the 
form 


0  =  6-  Acosw  —  ~  'T  (x-J--|£-)cos2co-r ... 


(3*0 


The  normals  to  the  pinch  and  housing  surfaces  in  the  same  approximation 
have  the  components 


n0  =  1 ,  /ta  =  —  sin  2w,  n9  =  O, 


(35) 


1  — 
A  . 


A* 


4p*(*  —  cos  2ito),  = 


—  — y  sw  to  +  2  y  sin  2«#  n9  =  0. 


(36) 


respectively. 

With  the  fulfullment  of  the  conditions  of  ideal  conductivity  of  the 
housing 


(37) 

(;t  is  the  discharge  time,  and  d  are  the  conductivity  and  thickness 
of  the  housing  wall)  the  normal  component  of  the  magnetic  field,: 
created  by  currents  in  the  plasma,  can  be  taken  a3  equal  to  zero  on 
the  housing.  Prom  thi3  condition,  A  and  6  are  determined.  The 

i  .  i 

influence  of  an  additional  magnetic  field  in  the  discharge  chamber 
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can  be  calculated  in  the  same  manner.  In  this  approximation  and  pro¬ 
vided  the  symmetry  IS  maintained,  the  additional  field  can  be  written 
as 

B  =  V<1>,  0=  Z?#gsin<i>  + 

+  B1^sin2a>=Btz  +  Blt^-,  &Q) 


where  B0  is  the  value  of  the  magnetic  field  component  normal  to  the 
housing  surface  in  the  center  of  its  section:  Bx  characterizes  the 
"barrel-shape  quality"  of  the  additional  field.  Disregarding  whether 
the  field  extends  to  the  housing  or  only  to  the  conductors  situated 
between  the  plasma  and  the  housing,  an  effective  boundary  condition 
will  be  the  equality  between  the  normal  component  of  field  B  outside 
the  plasma,  calculated  without  such  conductors,  and  the  value  of  the 
normal  component  of  the  additional  field  (38)  extended  to  the  housing 
i.e.  i 


(Bn)  =  B#sin  o>  ^  B,  —  sin  2<o. 


(59) 


The  magnetic  field  outside  of  the  plasma  is  determined  from  the  equa¬ 
tions 


B  =  VO;  } 

,  1  OO  1  <?»4)  __ 

7  0  a0  +  c*  da* ' 

__  sino)  OO  coso)  d<D  J 


Qda 


dQ 


On  the  pinch  surface  by  definition1 


(40) 


(Bn)  =>0.  (41) 

The  value,  of  the  co-component  Of  the  magnetic  field  on  the  pinch  3ur-  . 
face  we  will  write  as  '  ; 

I 

I 

-  1 

«  Ba (\  r~  A, cos u  r  -*V A*  cos 2<o). 

V  N  '  ;cr 


1  (42) 
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In  agreement  with  (21) ,  the  value  is 


8.*t  p 


f  ~  —  1 
i  2  i. 


(*3) 


The  precise  value  of  Aa  for  further  examination  is  immaterial. 

,  The  calculation  of  the  magnetic  field  B,  displacement  A,  and 
distortion  6  are  elementary;  thus,  below,  we  will  introduce  only  the 
calculation  result. 

1.  Magnetic  Field  Outside  Plasma.  In  the  examined  approximation 
the  solution  of  the  equations  in  (40) ,  on  fulfilling  boundary  condi¬ 
tions  in  (4l)  and  (42) ,  has  on  the  pinch  surface  the  following  form 


B0 


Qd  <i>  * 


where 

O  =  <I>ga>  -f  <Dj  sin  &>  +  sin  2«. 

Here  i 


=*  Bji\ 

m 

=  TR  B'‘a  [(Al  +  I)  (^+|)  + 

(46) 

“  y(A» -*•?)- + 

~r  "4o»"  (  A,  —  -j--r  + 

+  ^i‘CA*  +  '4 

1 

(*7) 

2.  Displacement 

of  f 

a  Pinch  in  a  Ideally  Conducting  Housing. 

The  displacement  A  in 

terms  of  order  ba/Ha) 

the  second  approximation  (with  respect  to 

! 

agrees  with  the  displacement  calculated  in 

the 
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3econd  approximation  and  is  determined  by  the  formula* 


.A 

6 


b 

'J.K 


(•-$)(*+*)]-*• 


Here  =  Ba  ^  At  is  determined  with  respect  to  (4j)  .  The 

relation  B0/B^  agreeing  with  the  usual  reference  frame  is  positive 
when  the  direction  of  the  additional  magnetic  field  B0  agrees  with 
the  direction  of  the  field  comprised  of  a  current  on  the  outside  of 

the  pinch.  . 

In  the  "Tokamuk"  device,  between  plasma  and  housing  there  is  a  1 
weak  conducting  chamber  made  of  thin  rust-proof  steel  through  which 
flows  the  current  equalizing  the  electrical  field  along  the  length  of 
the  discharge  chamber.  Let  bt  be  the  radiu3  of  the  chamber,  Ai  be  the 
displacement  of  the  center  of  the  chamber  cross  section  relative  to  , 
the  center  of  the  housing  section  in  a  direction  from  the  axis  of 

j 

symmetry.  Further  let  the  resistance  of  the  chamber  in  the  section 
change  with  respect  to  the  law  1  +  X  C03  o  (for  a  bellows  chamber  ; 
X  =  0,  for  smooth  chamber  X  ~  jr*-) ,  Then  the  additional  magnetic  field 
created  in  the  discharge  chamber  by  current  flowing  along  the 
chamber,  is  | 


■  r  i 

“V* 


(*9)  ! 


For  the  parameters  of  the  wTokamak-2"  device  £  14 J  (2R  =  125#  hi  =  20, 
b  =  25)  the  total  displacement  of  the  pinch  (in  centimeters)  with  a 


*  For  the  idealized  ca3e  of  surface  current  and  when  B0  =  0  the  j 
corresponding  calculation  of  displacement  in  the  first  approximation  ; 
ia3  done  by  V.  M.  Novikov  in  a  thesis  (1958)  completed  under  the,  i 

supervision  of  thi3  article* 3  author.  The  expression  for  A  where  , 
a  «  b  and  l  =  0  is  introduced  in  another  article  (11).  '  — 1 

0  1  . . . 9  ...  J 

°  •  ■*  *  f  s yj? 
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radius  equal  to  the  radius  of  the  diaphragm  (a  •»  10) ,  determinable 
(disregarding  Joints)  by  the  formulas  in  (48) ,  (49) ,  is 

A  =  2,5  +  4,2  + 1,9  A.  (50) 

The  additional  displacement  associated  with  the  currents  in  the  chamber 
cam  be  eliminated  by  selecting  appropriate  values  of  X  and  At.  When 
At  *  0  and  b  -  bt  «  bt  it  is  necessary  to  have  X  »  (this  condition 
is  fulfilled  when  the  chamber  is  smooth) .  For  the  bellows  chamber 
the  additional  displacement  can  be  reduced  to  zero  after  having  shifte< 
the  chamber  to  the  outside  housing  wall  to  a  distance  of  A  t  *  ^  (b  - 

-bt).* 

More  hazardous  is  the  shift  associated  with  a  pressure  increase 

while  heating  the  plasma;  this  shift  with  the  parameters  of  the 

,,Tokamak-2,,  reaches,  as  is  obvious  from  (50) ,  several  centimeters 

—  }1% 

when  the  pressure  changes  from  zero  to 

3.  Distorted  Shapes  of  a  Pinch  Section.  The  distortion  ampli¬ 
tude  6  is  determined  by  the  relation  * 


r— [(**)•+ 

2(A,+  *)l»±  +  (.— £)(A!+i.+ 

+v- 


(51) 


*  These  conditions  were  obtained  by  A.  I.  Morozov  and  L.  I. 
Solov*ev. 

**  We  notice  that  with  a  small  plasma  density  the  equilibrium 
conditions  can  be  influenced  by  the  centrifugal  force  of  the  current. 

for  its  calculation,  we  add  to  Ai  the  term  where  N  » 

Jt 

*  (when  v^  ■  const,  the  value  N  agrees  with  the  linear 

number  of  discharges  in  the  pinch) . 
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where  6  is  the  distortion,  caused  by  the  supplementary  scattered  fields 


<>*  )  lir«7v,nT- 

__  •••  a,  | 

2  J  K  th;  J- 


(52) 


The  negative  value  of  6  corresponds  to  elongation  of  the  pinch 
ection  along  the  axi3  of  symmetry  of  the  torus.  When  Bi  =  0  the 
expression  for  6  can  be  approximately  written  as 


6  2 A* 

—  rsj _ i _ 1 *  _ 

«  6*  * 

i  i 


Thus  it  is  obvious  that  with  a  limited  displacement  value  of  (A  <  b  -  a) 
the  comparative  distortion  is  always  very  small  (6/a  <  1/8) .  A  notice- 
able  distortion  of  the  section  is  possible  with  a  sufficiently  large 
constituent  of  the  quadrupole  component  of  the  additional  field  (Bt  *» 

~  B0) .  The  influence  of  this  field  when  B0  =  0  is  illustrated  in 

a  ( 

Fig.  2.  For  the  case  shown  in  Fig.  2  (Bi  >  0) ,  this  field  is  added 
to  the  current  field  on  the  lateral  sides  of  the  pinch  and  is  sub¬ 
tracted  from  it  on  the  upper  and  lower  part3  of  the  pinch. 

With  this  orientation  of  the  fields  the  pinch  will  be  stretched 
out  '’long  the  axis  of  symmetry.  Outside  of  the  plasma  the  general 
magnetic  field  change  associated  with  distortion  of  the  section  is 
determined  from  Formulas  («)-(4T)  and  has  the  form  (o>- component) 

bB„  =  Ba  A  (  I cos  2®.  (5*0 


4 .  F.ffect  of  Joints  on  Pinch  Displacement.  When  the 
a  joint  the  normal  component  of  the  magnetic  field  will  be 
tion  of  <p.  In  the  first  approximation 

i 

•  •  ;  2  ! 

(Bn)  =»  /  (<p)  sin  m. 

I  Vo  .  .  . 1  .  . 

;■  Sr  *£  Htvt* 


housing  has 
the  func- 

(55) 

i 
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I 


We  will  expand  f(<p)  in  a  Fourier  series  i 

■  i  .  : 

•  «) 

/(<p)  =  fl«+  2  V  cos m<p.  (56) 

The  harmonics  with  m  1  lead  to  the  appearance  of  heterogenity  (which 
we  will  not  calculate  here)  of  the  pinch  along  its  length,  the  zero 
harmonic 

2*  ; 

(57) 

,  i :  .t  t'  •'  iti.e  j 

leads  to  a  certain  displacement  of  the  column  as  a  whole.  For  a  zero': 

| 

harmonic,  (55)  has  the  form 

( Bn)  =  Ba  sinto,  (55a) 

i.e.,  it  agrees  with  the  condition  on  the  housing  (39)  in  the  presence 
of  the  additional  field  B0.  The  value  f(<p)  (in  the  absence  of  some 
kind  of  screen)  wa3  determined  with  logarithmic  accuracy  in  a  previous 
article  [10].  In  our  notations  it  is  written  asj 


<i>,  (b) 


/(< P)  = 


2 U  y\  —  yx  In— r- 

•  A 


0. 


^  .  (  -  <Pi  <  <P  <  ?P») 
(VXPI). 


The  width  of  the  Joint  is  taken  a3  equal  to  2h;  a  is  a  coefficient  of 

t 

the  order  of  unity  [10]j  <Pi  =  <&i(p)  sin  o>  is  the  correction  for  I. 

the  magnetic  "potential"  in  the  first  approximation  with  an  uninter¬ 
rupted  casing,  i.e.,  in  the  examined  case  4>i(p)  is  determined  by  ! 
Formula  (46) .  From  (57) »  (58)  we  will  find 


1 

0 


»Pi  «>) 
2 .iR 


► 


1  (59) 


1 


;  !CP 


-  wj? 


i 
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•  r 


Here  ^ : 


2  Jn 


ab  •  In  the  case  of  N  Joints  separated  from  one  another 


k  i 

by  a  distance  larger  than  b,  we  can  assume 

I 

J¥ 

U-V. 


21b 


0.6 


(60) 


Substitution  of  (59)  into  the  displacement  formula  in  (48)  gives 

+  (*+tX‘-«K*+'‘TH-)  +  (6!) 

+  [(2A1  +  .)-J-«]|1^r}. 

j 

Using  Conditions  (55a)  and  (57)  analogously,  we  can  refine  Expression 
(49)  for  the  field  created  by  the  current  in  chamber  Ji,  with  con¬ 
sideration  of  Joints  in  the  housing.  The  correction,  obviously,  has 

l  i 

the  order  n 


2 


5.  Pinch  in  a  Poorly  Conducting  Housing.  Having  condition  ! 

opposite  to  (57)#  i.e.,  in  the  case  where  the  thickness  of  the  skin-  | 
layer  is  larger  than  the  thickness  of  the  housing  walls,  the  asym¬ 
metrical  (with  respect  to  u>)  portion  of  the  magnetic  field  on  the 
housing  must  satisfy  the  condition: [10] : 

i 

(6j>)  ■ 

where  H®  is  the  field  outside  the  housing;  H^  is  the  field  between  the 
housing  and  the  plasma.  It  is  convenient  to  define  these  fields  in  ! 
terms  of  the  magnetic  current  by  Formula  (25) •  Considering  that 


E _ 

txcr  dt 

the  housing: 


from  (62)  we  will  obtain  the  following  condition  at 


Oq 


dQ 


4.*t  r/<T„ 

c*  "5 r 


(63) 


t  l 

»  1 


(with  p  we  designate  the  value  pe  =»  p1  at  the  housing)  . 


■.1 


hcr'b 
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<») 


4  6 

The  expressions  for  y  and  y  are  determined  by  the  formula  in 
(26)  (in  order  to  obtain  ye ,  we  take  a0  *  ai  *  0) 

+(a,+1)(1-£)]  “““}• 

— F  (.-'•  +  y)(,_>)]  "*  *}• 


(65) 


If  we  substitute  (64)  and  (65)  into  (63)  we  will  obtain  the  equation 
for  displacement  A 

-tin A-AJi-  - 

-«2bO"“+A--i)A-  •  (66) : 

where  A0  is  the  displacement  in  an  ideally  conducting  housing  (see  («)> 
where  B0  =  0: 


(6r) 


The  characteristic  time  of  holding  plasma  In  the  housing,  which  follows 
from  (66),  i3  determined  by  the  following  expression  (compare  with 

til)): 


_  A  M  -  ,  _  iY' 

*rt  ~  c*  ^  u  4  1  A*  2  )  * 


(68) 


In  conclusion  I  wish  to  thank  N.  A.  Yavlin3kiy  and  his  co-workers, 
especially  K.  A.  Razumova  for  the  interest  shown  in  these  calculations. 
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Fig.  1.  The  volume  element, 
from  whose  equilibrium  con¬ 
dition  is  determined  the  dis¬ 
tribution  of  the  magnetic  field 
on  the  pinch  surface. 


Fig.  2.  Configuration  of  the 
field  causing  distortion  of  the 
pinch  section. 
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